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, The percolation of Kandel, Ben-Av and Domany clusters for 2d fully frus- 

trated Ising model is extensively studied through numerical simulations. Crit- 
ical exponents, cluster distribution and fractal dimension of percolative cluster 
are given. 
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I. INTRODUCTION 



A 2d fully frustrated (FF) Ising model is a model with Ising spins ±1 where the interac- 
tions between nearest neighbor spins have modulus J > and sign ±1 (ferro / antiferromag- 
netic interactions) and where the signs are chosen in such a way that every plaquette (i.e. 
the elementary cell of square lattice) is frustrated, i.e. every plaquette has an odd number of 
-1 interactions so that the four spins of the plaquette cannot simultaneously satisfy all four 
interactions. In Fig. 1.1 we give an example of such a deterministic interaction configuration. 
In a plaquette of the FF model we can have only one or three satisfied interactions. The FF 
model has an analytical solution H and a critical temperature at T c = 0. 

Since single-spin dynamics for FF suffers critical slowing down, a fast cluster dynamics 
was introduced by Kandel, Ben-Av and Domany (KBD) in Ref. @. 

The KBD-clusters are defined by choosing stochastically on each plaquette of a checker- 
board partition of a square lattice one bond configuration between the three shown in Fig. 1.2. 

The probability of choice depends on spins configuration on the plaquette and it is a func- 
tion of temperature (correlated site-bond percolation ||). When there is only one satisfied 
interaction the zero-bond configuration is chosen with probability one. When three inter- 
actions are satisfied the zero-bond configuration is chosen with probability P = e~ iJ ^ kT ^ 
(where k is the Boltzmann constant and T the absolute temperature), the bond configuration 
with two parallel bonds on two satisfied interactions is chosen with probability P\ = 1 — Pq 
and the third bond configuration has zero probability. Two sites are in the same cluster if 
they are connected by bonds. For sake of simplicity from now on we choice J/(kT) = 1/T. 

Ref. || has stimulated several works j|,[5| that pay attention mainly to dynamics and 
to number of clusters and clusters sizes. In || numerical simulations on relatively large FF 
lattice sizes (number of sites N = 60 2 4- 120 2 ) supported the idea that the KBD-clusters 
represent spin correlated regions (like Coniglio-Klein clusters § do in Ising model) and 
consequently percolation temperature T p coincides with critical temperature T c , percolation 
exponents coincide with critical ones and KBD-clusters at T p are 2d self-avoiding walks 
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(SAW) at 9 point. 

In this paper we extensively study percolative features of KBD-clusters, considering very 
large lattice sizes (N = 100 2 -r-400 2 ), and give numerical results on critical exponents, cluster 
distribution and fractal dimension at percolation point. 

II. CRITICAL EXPONENTS AND PERCOLATION POINT 

We consider finite systems with increasing size (L = 100 -j- 400) with periodic boundary 
conditions. 

A cluster percolates when it connects two opposed system sides. For every size L there 
is a percolation temperature T P (L). With T p (without any argument) we mean the perco- 
lation temperature in the thermodynamic limit i.e. T P (L) — > T p for L — > oo. In this limit 
percolating clusters are present at T < T p but not at T > T p . In Fig. 2 we show typical 
clusters at several temperatures for a finite system with size L = 60. 

For every L we have studied the mean cluster size S = J2s s<2n s (where s is the cluster 
size, n s the number of cluster of size s per lattice site and the sum is extended over all finite 
clusters), the percolation probability P = 1 — J2 s sn s, the number of cluster N c = J2s n s, 
the number of bonds per lattice site A^, the mean size of the largest (percolating) cluster 
Si and the mean size of the second largest (percolating) cluster Sjj. These quantities are 
shown in Fig.3 for L = 100 -r 400. Let's note that for T — > the bonds cover 50% of lattice 
interactions (that is the random-bond percolation threshold on square lattice), Sjj goes to a 
finite value (like predicted by KBD || and already verified in Ref. 0) and occupies almost 
35% of the lattice, and that Sj occupies almost 65% of the lattice. At T = only two 
clusters survive, as shown in Fig. 2 e) and f). 

Now we will give numerical estimates of critical exponents that characterize the KBD- 
cluster percolation. 

We know || that in the thermodynamic limit the mean cluster size diverges for T — > T p , 
the percolation probability goes to zero in the limit T — > T~ and the number of cluster goes 



3 



to zero for T — > T p . 

We assume that near T p the connectivity length £ (i.e. the typical linear cluster size) 
diverges like £ ~ \e~ 2 l T — e~ 2 l Tv \~ u , the mean cluster size diverges like S ~ |e~ 2//T — e~ 2 / Tp |~ 7 , 
the percolation probability goes to zero like P ~ |e~ 2 / T — e~ 2//Tp I' 3 and the number of cluster 
goes to zero like iV c ~ |e~ 2 / T — e _2//Tp | 2_a . The last relations are definitions of critical 
exponents a, (3, 7 and za 

By standard finite-size scaling considerations |§ we can make the ansatz 

S~L^f s (\e- 2 / T -e- 2 ^\L^) (1) 

P ~ L-^/p(|e- 2/T - e-^lL 1 /-) (2) 

and 

iV c ~ L^/"f Nc (\e- 2/T - e-^'lL 1 '") (3) 

where fs{x), fp{x) and fN c ( x ) are universal functions, i.e. independent by L. 

Via data-collapse (see Fig.4) we estimate the parameters e~ 2 ^ Tp = 0.0000, a = 0.1, 
(3 = 0.00, 7 = 2.00, and v = 1.00 with error of one unit in the last given digit . Therefore 
the scaling relation a + 2(3 + 7 = 2 and the hyperscaling relation 1 — a = vd are satisfied 
with good approximation. 

In Tab.l we give numerical estimates of T p (L). The data are obtained taking for L = 
100,200,300,400 the values of T p (L) at which the S data in a log-log plot vs. \e~ 2 / T - 
e - 2 / T p( L )| follow two parallel straight lines (one above and one below T p (L)) with slopes in 
good agreement with 7 = 2 and then best-fitting these values as e~ 2 ^ Tp ^ ~ 1/L. 

III. FRACTAL DIMENSION AND CLUSTER DISTRIBUTION 

Let's now consider the fractal dimension D of the percolating cluster. From the scaling 
invariance hypothesis || we know that P ~ £, D ~ d , then we obtain D = d—(3/v (hyperscaling). 
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In present case we have (3 = 0, then D = d = 2. The same result is obtained from the scaling 
relation (3 + 7 = Dv. 

This is confirmed by the analysis of cluster distribution (see Fig. 5). The scaling invariance 
hypothesis |||9| gives for T — > T p and s — > 00 

n, = S -Vn s (|e- 2 / T -e- 2 ^| S CT ) (4) 

with r = 1 + d/D, o = 1/{vD) and /n s (x) universal function. 

From data-collapse for n s near T p (see Fig. 6. a) we obtain numerical estimates of param- 
eters. The data in Fig. 6. a are chosen in such a way that the quantity (e~ 2 / T — e~ 2 ^ Tp )L 1 ^ u 
(with e~ 2 l T v = and v = 1) is a constant with T ~ T P (L) for every considered L. The 
results are r = 2.00 and o = 0.50 (with error of one unit in the last digit), that, with the 
definitions of r and a, give D = 2 and v = 1. On the other hand these values of r and cr 
satisfies the relation a{2 — a) = r — 1, a/3 = r — 2, (77 = 3 — r. 0] 

From Fig.6.a we see that the universal function f ns (x) is a bell-shaped curve for 
T ~ T p (L). For temperatures slightly below T p (L) (Fig.6.b) f ns (x) is shifted, while for 
temperatures slightly above T p (L) (Fig.6.c) f ns (x) changes dramatically its shape. 

Away from T P (L) we know || that is valid the relation 

logn s ~ -s c (5) 

for s — >• 00, with ( = 1 above T p (L) and C = i — 1/d = 1/2 below T p (L). This relation is 
confirmed with reasonable approximation by our numerical simulations, as shown in Fig. 7. 
Let's note that, while the exponent ( = 1 above T P (L) is good for a wide range of s 
(s = 2000 -7- 8000 for L = 100), the exponent ( = 1/2 below T P (L) is good for a smaller s 
range (s = 2000 -r- 4400 for L = 100) since finite-size effect become more important below 
T p (L). The smaller T, the smaller s range is. 

A direct way to estimate the fractal dimension D is given through its definition 

s ~ R D (6) 
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for T = T p , with R radius of gyration of the cluster of size s. We know that cluster 
dimension deviates from D away from T p , becoming the Euclidean dimension d below T p and 
a value smaller than D above T p . This is true because eq.(^) is valid within the connectivity 
length £ for all temperatures, but £ goes to zero away from T p . Unfortunately data about 
this relation are difficult to analyze. Indeed near T p (L) for every finite system with L < 120 
it seems that D is almost 7/4 = 1.75 (the fractal dimension of a SAW at 9 point), but for 
larger L (see Fig. 8 and Tab.l) the fractal dimension D grows slowly to the asymptotic value 
2. 



IV. CONCLUSIONS 

We have numerically investigated the KBD-cluster percolation problem in 2d FF Ising 
model. From our simulation we found that, within numerical errors, this correlated site-bond 
percolation satisfies scaling and hyperscaling relations and have, in the thermodynamical 
limit, a percolation temperature T p = and the exponents a — 0, (3 — 0, 7 = 2, v — 1, 
T = 2, a = 1/2, C(T > T P (L)) = 1, ({T < T P {L)) = 1/2. Moreover at T p clusters are 
compact (fractal dimension D = 2). Therefore now we can correct the conclusion of Ref. 
H and say that, since T p = T c and v is equal to spin correlation exponent [10|, the site 



connectivity length £ goes like the spin correlation length diverging at zero temperature. 
Although the exponent 7 is different the coincidence between £ and correlation length is 
enough to give an efficient Monte Carlo cluster dynamics. 
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TABLES 



TABLE I. Numerical estimates of T P (L) and D(L) for L = 60-^400. The way used to evaluate 
T P (L) and D(L) give us confidence only on digit not in parentheses. 



L 


60 


80 


100 


120 


200 


300 


400 


T P (L) 


0.51(7) 


0.48(1) 


0.45(6) 


0.43(7) 


0.39(2) 


0.36(1) 


0.34(2) 


D{L) 


1.7(2) 


1.7(5) 


1.7(7) 


1.7(9) 


1.8(2) 


1.8(5) 


1.8(6) 
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FIGURES 

1) 



2) 



a 



b 



c 



FIG. 1. 1) Example of 2d FF lattice: the spin are on the vertices; solid lines represent 
ferromagnetic interactions (+ J) and dashed lines antiferromagnetic interactions ( — J). 2) Plaquette 
bond configurations: a) zero bond; b) two parallel vertical bonds; c) two parallel horizontal bonds. 
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a 



b 




FIG. 2. Typical KBD-clusters on a FF lattice with size L = 60 with periodic boundary 
conditions: a) at T = 1; b) at T = 0.65; c) at T = 0.53 slightly above T P (L) ~ 0.52; d) at 
T = 0.52 ~ T p (L); e) at T ~ (largest cluster); f ) at T ~ (the second cluster). 
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FIG. 3. Mean cluster size S, percolation probability P, number of cluster N c , number of bonds 
per lattice site N^, mean size of largest cluster Si and mean size of second largest cluster Su vs. 
temperature T for lattice sizes L = 100, 200, 300, 400. The error bars (often included in symbols) 
are the statistical errors. 
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FIG. 4. Scaling for S 1 , P and iV c following assumptions @), ((H) and (||) for data of systems with 
sizes L = 100, 200, 300, 400. The parameters e~ 2 ^ = 0.0000±0.0001, a = 0.1±0.1, = O.OOiO.01, 
7 = 2.00 ± 0.01 and v = 1.00 ± 0.01 are such that the data for different sizes L collapse on single 
curves (one for each graph). These curves are, respectively, the universal functions /s, fp and fjy c . 
The errors are estimated observing the range of parameters within which the data points collapse 
roughly well. 
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FIG. 5. Cluster distribution n s vs. cluster size s for a system with size L = 100 at several 
temperatures T. Every bin is large 400 unities in cluster size. The percolation temperature for this 
system size is T P (L = 100) ~ 0.46 then for all T > 0.47, even if the highest bin is not empty, there 
are only no percolating clusters. Let's note that for all T above T P (L) it is n s ~ e~ s and that the 
distribution become symmetric for T — > 0. 



13 



o 
o 



□ 



o 
o 

C\| 



o 
o 



o 
o 




s u 



CO 



d 



d 



d 



o 




FIG. 6. a) Scaling for n s following assumption (0) with parameter r = 2.00 ± 0.01 and 
a = 0.50 ± 0.01 for data of systems with sizes L = 100, 200, 300, 400. Each set of data is cho- 
sen at a temperature near the corresponding T p {L). As consequence for each temperature the 
quantity (er 2 l kT — e~ 2 / Tp )L 1 / u with T p = and v = 1 is equal to 2.084. Every point in the graph 
is an average over 500 consecutive values of s for L = 100, 2000 for L = 200, 4500 for L = 300, 
8000 for L = 400. b) As in part a) but for (e" 2 / T - er' 2 l T p)L 1 l v = 0.823 (below T p (L)). c) As in 
part a) but for (e" 2 / T - e^^^L 1 ^ = 3.567 (above T P (L)). 



14 




1 1.2 0.7 0.8 0.9 1 

log (s/400) log (s/400) 



FIG. 7. Fit of log(-logn s L 2 ) vs. log(s/400) for a system with L = 100 and T P (L) ~ 0.46: a) 
at T ~ 0.47 > T P (L) the slope is C — 1; b) at T ~ 0.45 < T P (L) the slope is C — 1/2. 
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FIG. 8. Estimate of D from definition (§) vs. l/L (see Tab.l): at T = 0.517 for L = 60, at 
T = 0.481 for L = 80, at T = 0.450 for L = 100, at T = 0.437 for L = 120, at T = 0.389 for 
L = 200, at T = 0.360 for L = 300 and at T = 0.343 for L = 400. The error bars are probably 
underestimated. The arrow heads for the asymptotic value of D. 
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